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Abstract 

We have previously developed a chemical conversion model of the carbothermal process- 
ing of lunar regolith using methane to predict the rate of production of carbon monoxide. 
In this carbothermal process, gaseous methane is pyrolyzed as it flows over the hot surface 
of a molten zone of lunar regolith and is converted to carbon and hydrogen. Hydrogen is 
carried away by the exiting stream of gases and carbon is deposited on the melt surface. 
The deposited carbon mixes with the melt and reacts with the metal oxides in it to produce 
carbon monoxide that bubbles out of the melt. In our model, we assume that the flux of 
carbon deposited is equal to the product of the surface reaction rate constant 7 and the 
concentration of methane adjacent to the melt surface. Similarly, the rate of consumption 
of carbon per unit volume in the melt is equal to the product of the melt reaction rate 
constant k and the concentrations of carbon and metal oxide in the melt. 

In this paper, we describe our effort to determine 7 and k by comparison of the predictions 
from our model with test data obtained by ORBITEC (Orbital Technologies Corporation). 
The concentration of methane adjacent to the melt surface is a necessary input to the 
model. It is inferred from the test data by a mass balance of methane, adopting the usual 
assumptions of the continuously-stirred-tank-reactor model, whereby the average concen- 
tration of a given gaseous species equals its exit concentration. The reaction rates 7 and k 
have been determined by a non-linear least-squares fit to the test data for the production 
of carbon monoxide and the fraction of the incoming methane that is converted. 

The comparison of test data with our model predictions using the determined chemical 
kinetic rate constants provides a consistent interpretation of the process over the full range 
of temperatures, pressures, and methane flow rates used in the tests, thereby increasing 
our confidence to use the model for scale-up purposes. 


Nomenclature 


a Stoichiometric coefficient 

A Melt surface area, m 2 

C Concentration, moles/m 3 

C m Concentration of methane in the gas phase adjacent to the melt surface, moles/m 3 
/ Fraction of incoming methane converted 

k Melt reaction rate constant, m 3 / (moles ■ s) 
m Molar flow rate, moles/s 

t Scaled (dimensionless) time, t = Yfkfo 

t* Time, s 

V Melt volume, m 3 
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X 

X 

y 

z 

Greek 
a 
a * 

P 

P* 

Pc 

7 

A 

P 


Mole fraction 
Number of moles 

Quantity defined as y = aJ y^ 2 m = a ^yjp in , m 3 / (moles ■ s) 
conversion fraction of metal oxide, z = a y%° 


Scaled concentration of carbon in the melt, a = 
Concentration of carbon in the melt, moles/m 3 
Scaled concentration of metal-oxide in the melt, p = ^ 
Concentration of metal-oxide in the melt, moles/m 3 
Initial concentration of metal-oxide in the melt, moles/m 3 
Surface reaction rate constant for methane, m/s 
Dimensionless parameter, A = 




Molar density of the gas phase, moles/m 3 


Subscript 

co carbon monoxide 
g gas phase 

h hydrogen 

in into the reactor 
m methane 
out out of the reactor 


I. Introduction 

A key element in NASA’s plans for exploration of the Moon involves the production of oxygen for life 
support as well as use as a propellant. Sanders et al . 1 describe NASA’s current plans for In-Situ Resource 
Utilization (ISRU). Various processes are being considered for the extraction of oxygen that is available in 
the lunar regolith as oxides of various materials. These processes include reduction using hydrogen, molten 
salt electrolysis, and carbothermal processing. 

The carbothermal processing of lunar regolith using methane as the source for carbon was first performed 
by Rosenberg and co-workers at Aerojet in the 1960s (Rosenberg et al .). 2 Currently it is being developed 
further for lunar exploration by Orbital Technologies Corporation (ORBITEC ). 3 We have previously devel- 
oped a chemical conversion model of the carbothermal processing of lunar regolith using methane to predict 
the rate of production of carbon monoxide (Balasubramaniam et al .) 4 In this process, gaseous methane is 
pyrolyzed above the hot surface of a molten zone of lunar regolith, depositing carbon on the surface and 
generating hydrogen which is carried away by the exiting stream of gas. The deposited carbon mixes with 
the melt and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. 
Oxygen is ultimately produced by processing the carbon monoxide in reactors downstream. In our model, 
we assume that the flux of carbon deposited (moles/ (m 2 s)) is equal to the product of the surface reaction 
rate constant 7 (m/s) and the concentration of methane (moles/m 3 ) adjacent to the melt surface. The rate 
of consumption of carbon per unit volume in the melt is equal to the product of the melt reaction rate 
constant k (m 3 /(moles -s)) and the concentrations of carbon and metal oxide in the melt (moles/m 3 ). 

In this paper we describe our effort to determine 7 and k by comparison of the predictions from our model 
with data from ORBITEC’s Test Series- 1. 5 The data from Test Series -1 is reproduced in Table 1 . This series 
of tests was conducted using the lunar regolith simulant JSC -1 A. A 200W laser was used to melt the material 
and maintain it at the chosen operating temperature. This test series was a part of ORBITEC’s Statistical 
Experiment Design to investigate the main effects of processing pressure, temperature, and methane flow 
rate on the carbothermal reduction process, as well as the interactions that occur among these test variables. 
The chemical processing time in these tests was 80 min. 
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II. Model 


A. Background 

Figure 1 shows a schematic of the carbothermal process. A model of the processing of lunar regolith by this 
method is described in Ref [4] . It is assumed that carbon deposition occurs on the hot molten surface via a 


Heat flux 



Figure 1. Sketch of the carbothermal process 

surface reaction. This surface reaction is regarded as first order with respect to a gaseous carbon- containing 
species and produces solid carbon and gaseous hydrogen. This carbon-containing species need not necessarily 
be methane; it can be any of several carbon- containing species which may be produced by the gas-phase 
pyrolysis of methane. Since all these carbon-containing species are ultimately produced from methane, the 
model lumps all possible pathways to deposit carbon on the surface into a single process. This deposition 
process is described by a first-order surface reaction of methane. Thus, the flux of carbon deposited on the 
surface of the melt is assumed to be proportional to the concentration of methane in the gas phase above 
the melt. The proportionality constant needs to be determined by comparison of the model predictions with 
the test data. The concentration of methane above the melt is assumed constant by its continuous flow in 
the gas phase. The deposited carbon and the metal oxide in the molten regolith are assumed to undergo 
a single-step, second-order chemical reaction (first-order with respect to both carbon and the metal oxide). 
This reaction rate constant also needs to be determined by comparison of the model predictions with the 
test data. Our current model does not account for spatial gradients of the concentration of carbon or of 
the metal oxide within the melt; the zone is regarded as well-mixed. This would be the case when the zone 
is sufficiently small that the diffusion of carbon within the melt is rapid, or more likely, when there is flow 
within the melt caused by the evolution of bubbles, natural convection, and thermo capillary forces. In some 
of the tests in Ref [5], a crust or cap of carbon forms on the surface of the melt (see Table 1). The formation 
of a carbon cap indicates that the deposited carbon did not get mixed in the melt, or more likely, did not go 
into solution in the melt. We do not consider the formation of a carbon cap in our model. We assume that 
the molten zone is established and maintained by heating the regolith by suitable means. The heat transfer 
aspects of the problem are not analyzed here. The processes of chemical conversion and heat transfer within 
the molten regolith are assumed to be uncoupled. 

B. Reaction in the melt 

The chemical conversion model in Ref[4] assumes that the concentration of methane in the gas phase is 
maintained constant. Carbon is deposited on the melt by a first-order surface reaction. Within the melt the 
metal oxide and carbon are assumed to be spatially uniform, and undergo a single-step reaction that is of 
first order in each species. The mass balance of metal oxide and carbon were described as follows. 

^ = ^-[.• 1 3 * 

dt* V H 
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( 2 ) 


dfi* 

dt* 


= -a k a* (3* 


where V = melt volume (m 3 ), A = melt surface area (m 2 ), a = stoichiometric coefficient of the oxide term 
in its reaction with carbon (« 0.5 for Si02), and C m = concentration of methane in the gas phase adjacent 
to the melt surface (moles/m 3 ). The initial conditions are that at t* = 0, there is no carbon present in the 
melt and concentration of the metal oxide is /?o- A closed form solution can be obtained and is given below. 


0 = 


a = p — itAi 


2y/\ Exp ( t — ^i 2 ) 


2v^A + Exp ( A-) y/27r 


Erf (^) + Erf ( 


V2A 1 )] 


(3) 

(4) 


where a = and /3 = denote dimensionless concentrations of carbon and metal oxide, respectively, 
and t = 1 j* k denotes scaled time. A = a k y^T is a dimensionless parameter that is a measure of the 
relative strength of the surface reaction rate of methane to the reaction rate between carbon and metal oxide 
in the melt. 


C. Gas phase above the melt 

In order to apply the melt reaction model described above, the gas-phase concentration of methane C m 
must be known. If the gas phase is flooded with pure methane by maintaining a high flow rate so that the 
mole fraction of methane is nearly unity everywhere, then can be determined from the conditions of 
the incoming stream of methane and the temperature of the melt surface. If a significant fraction of the 
incoming methane is consumed by the surface reaction to deposit carbon on the melt, then C m must be 
determined by considering the transport of various species in the gas phase. 

We consider the mass balance of methane, hydrogen and carbon monoxide in the gas phase above the 
melt. (We neglect the production of carbon dioxide in the melt. The experiments at ORBITEC indeed 
show that the amount of carbon dioxide detected is only a few percent of the amount of carbon monox- 
ide produced.) Let subscripts m, h and co denote methane, hydrogen and carbon monoxide, respectively, 
and subscripts in and out denote quantities at the inlet and exit of the reactor. Neglecting diffusion, the 
gas-phase molar balance for these species can be written as 


(^m)m (jdjnjout — H“ Vg ^.* 

(5) 

-{m h )out = —27 AC m + V 9~^r 

( 6 ) 

^ \ Vfi ]k~ aiJ ,dC co 

(mco)out — a + Vg 

(7) 


where V g is the volume of the gas phase and m is the molar flow rate of a given species. The term involving 
7 in Eq(5) represents the consumption of methane by the chemical reaction at the melt surface. In writing 
the mass balance for hydrogen, we have made use of the fact that the number of moles of hydrogen produced 
is twice the number of moles of methane pyrolyzed, which follows from stoichiometry. The production term 
in Eq (7) represents the carbon monoxide produced by the reaction in the melt that bubbles out of the melt 
and is released into the gas phase. 


We assume that the temperature and pressure are constant in the gas phase, and the gases obey the ideal-gas 
law. Under these conditions, the total molar density of the gases (moles/m 3 ) p = C m -\-Ch-\-C co is a constant. 
The addition of the three equations given above yields 


(fh-m)i 




out 


-7 AC„ 


vplk 


a{3 


(8) 


where J out “I - (di*h ) out “I - ( 'dico')out‘ 
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We further assume that: 

(i) the reactor can be modeled as a continuously-stirred-tank-reactor (CSTR) which implies that the 
average concentration of a given species in the reactor equals its exit concentration. Thus, = 
(Cm) out = %mP where is the mole fraction of methane. 

(ii) a quasi-steady state holds, and the unsteady (accumulation) term in the mass balance equation for 
methane can be neglected. 

(iii) a fraction / (taken to be a constant) of the incoming methane is pyrolyzed. 

The mass balance for methane and the overall gas-phase mass balance can then be written as follows. 


a 


f {jh*rn)i 
f) (^h-m)in _ 


7 AC m = ryApXn 


(^m)in "h 7 ApX r 


. + Y&afi 


(9) 


( 10 ) 


We note that when the flow rate of methane is large compared to the rate at which it is consumed by the 
reaction on the melt surface and the rate of production of carbon monoxide (i.e., ^ C { n - , 7 - ^ — <C 1), Eqs 
(9) and ( 10 ) show that / <1 and x m — ► 1 . This represents the case where the gas phase above the melt is 
flooded by methane, as discussed earlier. 


D. Determination of k and 7 

It can be shown that the amount of carbon monoxide produced, Xqo (moles) is given by 


Xco 


V(3 0 


a -m) 


Substituting the expression for f3 from Eq (4), this equation may be written as 
Vg Exp (k/3 0 t* - |yk/3pt* 2 ) 


2VT + Ex P (£) [Erf (yj) + Erf (^(yA)t* - 1)) 


= 1-z 


( 11 ) 


(12) 


where y = aJ y^ rn = (m 3 / (moles -s)), z = conversion fraction of metal oxide = a yp° • We have 

substituted 7 ^ 4 C m = /(m m )^ n (from Eq (9)) in the definition of y to eliminate 7 in the expression for y. 
Note that y/k is a dimensionless quantity (y/k = A), y is proportional to the rate of carbon deposited per 
unit volume of the melt. 


The reaction rate constant k is obtained from the above expression by a non-linear least-squares fit. From 
the perspective of the fitting procedure, y is the variable, 1 — z is the data and k is the parameter to be 
obtained. 


From Eqs (3), (9) and (10), the following expression can be obtained for 7 . 

f(l + f + f 0 t)(m m ) in kVfil ( f \ 

7 w + 

Once k is determined, 7 can then be obtained from a least-squares fit to this expression. 


(13) 
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III. Results 


The Test Series-1 data 5 is tabulated in Table 1. The area of the melt-gas surface is necessary in analyzing 
the data. We have assumed that the processed regolith has a hemispherical shape to determine this area. 
We have also assumed that JSC-1A has a maximum oxygen yield of 25% in calculating /3q. This is consistent 
with the extraction of all the available oxygen from Si0 2 that JSC-1 A contains, which is 46.2% by weight. 6 
For another lunar simulant JSC-1, it has been estimated that the maximum oxygen yield is approximately 
28% for reduction of Si0 2 , FeO and Fe 2 0 3 to Si and Fe, and Ti0 2 to Ti 2 0 3 . 7 The initial concentration of 
the metal oxide in the simulant is assumed to be /3q. 

Table 1 . Summary of Results from Test Series #1 



Experimental Variables 

Results 

Additional Observations 

Test 

Number 

Target 

Processing 

Temperature 

(C) 

Processing 

Pressure 

(psia) 

Methane 
Flow Rate 
(seem) 

Oxygen 

Produced 

(g) 

Oxygen 

Yield 

(%) 

Carbon 

Recovered 

(%) 

Average 

Processing 

Temperature* 

(C) 

Maximum 

Reflector 

Temperature 

(C) 

Regolith 

Processed 

(g) 

Amount of 
Methane 
Pyrolyzed 

(%) 

Time to 
Form Carbon 
Cap (min.) 

1 

1850 

9 

6.05 

0.340 

12.6 

79.1 

1846 

609 

2.702 

78.1 

no cap 

2 

1650 

9 

11.56 

0.302 

13.3 

62.0 

1735 

597 

2.268 

65.8 

55.2 

3 

1850 

15 

11.56 

0.402 

10.5 

73.3 

1890 

578 

3.845 

74.5 

58.9 

4 

1850 

3 

17.18 

0.160 

2.0 

75.8 

1911 

601 

8.119 

59.6 

no cap 

5 

1850 

15 

11.56 

0.342 

9.5 

67.7 

1913 

525 

3.587 

71.6 

60.8 

6 

1750 

9 

17.18 

0.301 

11.0 

52.6 

1713 

636 

2.740 

56.1 

32.1 

7 

1650 

3 

6.05 

0.219 

8.7 

76.7 

1722 

596 

2.516 

56.7 

no cap 

8 

1650 

15 

17.18 

0.295 

15.6 

50.1 

1702 

612 

1.889 

54.9 

16.6 

9 

1750 

3 

11.56 

0.325 

7.6 

69.8 

1713 

613 

4.298 

59.1 

no cap 

10 

1750 

15 

6.05 

0.289 

8.2 

83.3 

1666** 

597 

3.514 

77.7 

no cap 

11 

1650 

3 

6.05 

0.235 

6.1 

79.9 

1640 

638 

3.860 

56.9 

no cap 


* = Average processing temperature is calculated before carbon cap formation. 


= The pyrometer may have had an obscured view of the melt in this test, so the temperatures may not be accurate. 


We have determined that k = 4.512 x 10~ 8 m 3 /(moles-s) and 7 = 0.0404 m/s are the best-fit values to 
the data in Table 1. The Test Series-1 data in Table 1 is plotted in the variables used in the statistical fit 
(1 — z versus y ) in Figure 2 and the best fit curve for k is shown. The conversion fraction z is the fraction 
of the initial metal oxide present in the simulant that has reacted. It is readily obtained as the ratio of 
the actual oxygen yield in each test to the maximum possible oxygen yield. The data in Figure 2 and all 
subsequent figures shown are color-coded as follows: green symbols are used to represent data in Table 1 
where no carbon cap is formed (test number 1, 4, 7, 9, 10, 11), with a darker green used for test number 4; 
blue symbols represent data where the carbon cap formation is relatively slow (test number 2, 3, 5) and red 
symbols are used when the carbon cap formation is relatively fast (test number 6, 8). 



Figure 2. Data and best-fit curve for the reaction rate from Test Series-1. The best-fit value is k = 4.512 x 

10 -8 m 3 /(moles ■ s). 


The comparison between the predicted and the measured conversion fraction of the metal oxides in the 
simulant is shown in Figure 3. The data together with the best-fit curve for 7 is shown in Figure 4. It is 
evident from these figures that the theoretical model does capture the trend in the data. There is some 
scatter in the data when compared to the predictions. Test number 4 in the plot for 7 (Figure 4) appears 
to be an outlier. Also, in Figures 2 and 3, the data for test number 4 has a trend that is different from the 
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Figure 3. 


Comparison of measured and predicted conversion fraction from Test Series- 1 . 


rest of the no-carbon-cap data. Indeed, the test data in Table 1 shows that for test number 4, the amount 
of regolith processed is high and the oxygen yield is low relative to the other tests. For these reasons, a dark 
green symbol has been used for test number 4 in all the plots. 


Y 


0.16 
0.14 
0 . 12 
0.1 
0.08 
0.06 
0.04 
0.02 


10 


Test Number 


Figure 4. Data and best fit curve for the pyrolysis constant from Test Series- 1. The best-fit value is 7 = 0.0404 

m/s. 


There appears to be a trend in Figures 2 and 3 with respect to carbon cap formation. A carbon cap does not 
occur when the conversion fraction of the metal oxides is low. When a carbon cap is formed, the conversion 
fraction is high, and is typically above 0.4. The trend is reasonable because, when the conversion fraction of 
the metal oxide is high, the amount of unreacted metal oxide in the melt is low, and therefore, the reaction 
rate in the melt is relatively slow. Thus, the amount of unreacted carbon in the melt would be high, and 
there is increased resistance for additional carbon that is deposited on the surface to go into solution in the 
melt. One might anticipate that there is a tendency for the carbon deposited on the surface to agglomerate 
and form a layer or cap on the surface. 

The analysis of the data that has been performed assumes k and 7 to be constants, independent of temper- 
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ature and pressure. In the Test Series 1 data, the processing temperature varies from 1640°C to 1913°C, 
and the pressure varies from 3 psi to 15 psi. These variations are accounted for in the data analysis by 
their effect on the gas density. One might speculate that there may be additional temperature and pressure 
dependencies that are not captured in the model if k and 7 were to depend on temperature and pressure. 
It is an open question whether the data scatter in the figures might be due to these temperature and pres- 
sure dependencies. Since the data in Table 1 is scarce, further test data is likely necessary to explore such 
dependencies. These features need to be addressed in the future. 

IV. Conclusion 

The comparison of ORBITEC’s Test Series 1 data with our model of the pyrolysis of methane and reaction 
of the deposited carbon with the metal oxide in the melt yields a surface reaction rate of 7 = 0.0404 m/s and 
a melt reaction rate of k = 4.512 x 10 ~ 8 m 3 /(moles • s). The data encompasses processing temperature from 
1640° C to 1913° C, processing pressure from 3 psi to 15 psi, methane flow rate from 6.05 seem to 17.18 seem, 
and a processing time of 80 min. With the best-fit values of 7 and /c, the model predictions for the amount 
of metal oxide that is converted compares favorably with the experimentally measured conversions over the 
entire range of test conditions. Thus the model can be used with a good degree of confidence to predict the 
behavior of larger scale carbothermal reactors. One aspect that needs to be understood is the scatter in the 
data when compared with the model predictions, and whether temperature and pressure dependencies of 7 
and k can explain or reduce the scatter. 
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